COMPUTERS 


Electro/|/Scientific Industries 


7524 S.W. MACADAM AVE. + PORTLAND 19, OREGON JANUARY 1961 


SAMPLED-DATA SYSTEM DESIGN WITH A POTENTIAL-PLANE 
ANALOG OF THE LOG Z-PLANE 
by 


James C. Looney 


F(z) 


e*(t) 
SAMPLED 


+ 
SAMPLING 
INTERVAL 


ee enn wer ie ee ere eee =— ss 


OUTPUT 


€ 
, CONTINUOUS 


ERROR —— 
e(t) COMPENSATION DATA CONTINUOUS 


NETWORK RECONSTRUCTION SYSTEM 


INPUT 


A SAMPLED-DATA CONTROL SYSTEM 


© Copyright 1961 by Electro Scientific Industries, Inc. el elsii 
LITHO IN U.S.A. 


SAMPLED-DATA SYSTEM DESIGN WITH A POTENTIAL-PLANE 
ANALOG OF THE LOG Z-PLANE 
INTRODUCTION 

Extensive use of digital computers in the past decade has established a 
definite trend toward control systems that operate with variables which are in 
the form of a sequence of numbers. These variables are usually quantized in 
amplitude and available only at specific, equally spaced instants of time. In 
comparison, the continuous system has variables that are continuous functions 
of time. In the systems using discrete variables, the data appear as sequences 
of numbers or pulses and these data are valid only at the sampling instants. This 
type of system is usually called a sampled-data system. 

Several methods of analyzing and synthesizing sampled-data systems have 
been discussed in the literature. One popular and effective method uses the Z- 
transform and root-locus: plots of the closed loop system on a geometric analog 
of the z-plane. This paper shows how sampled-data problems can be solved with 
a potential-plane analog of the log z-plane. This method combines the advantages 
of fast solutions from a potential-plane analog and easy visualization of system 


changes on the log z-plane. 
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THEORY 

One of the fundamental decisions to be made concerning a sampled-data sys- 
tem is what sampling period or frequency is to be used. Shannon has found that 
the theoretical minimum sampling frequency that can be used without loss of in- 
formation is equal to twice the value of the highest frequency of interest in the data 
to be sampled. However, the minimum sampling frequency is based upon the use of 
an unattainable ideal, low-pass filter to reconstruct the sampled data, so practical 
systems use a sampling frequency which is usually much higher than the theoretical 
minimum. 

Figure 1 shows a block diagram of a typical sampled-data feedback control 
system. The effect of the sampler is to convert a continuous function of time é (t) 
into a sequence of pulses e*(t) whose amplitudes are the same as e(t) at the samp- 


ling instants. 
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FIGURE | A SAMPLED-DATA CONTROL SYSTEM 


The relationship between any continuous function f(t) and the corresponding 
sampled function f*(t) has been derived by several authors (1,3,4) and only the re- 
sults will be repeated here. If the assumption is made that the width of the sampling 
pulse is small compared to the sampling interval T, the mathematics can be simpli- 
fied by the use of the impulse function S(t). 


For this case, 


e*(t) = £(1) Y 8(t=nT) (1) 
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ae 
Since the impulse is very narrow, the actual value of f(t) is f(nT), the value of the 


function at the sampling instant. Thus Equation (1) can be written as 
% ice) 
f*(t) = > f(nT) 8(t—nT) (2) 
n=O 
If the Laplace transform is taken, term by term, the result is 


F*(s) = 2 F(nT) gous (3) 


It has been shown that an equivalent expression for F*(s) can be obtained with the 


aid of a Fourier series for 8(t—nT) which is 


+00 
-  F¥(s) = 2 Y F(s+njws) (4) 
n=-0 
where 
ee 
a 


This expression shows that the effect of combining a sampler with F(s) results in 
a function F*(s) which is periodic in Jws. Figure 2 shows the s-plane plots of an 


F(s) and the corresponding F*(s). 
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FIGURE 2 THE s-PLANE PLOTS OF A SYSTEM RESPONSE; (a) WITHOUT A 
SAMPLER, (b) WITH A SAMPLER 
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It is easily seen that F*(s) consists of an infinite number of poles and zeros 
which come from the repeated horizontal strips of width jwWs. In practical prob- 
lems, F™(s) can be evaluated by ignoring the poles and zeros that are far removed 
from the ¢ axis. However, if F(s) contains many poles and zeros, the evaluation 
of F*(s) will be impractical with analog methods. 

The evaluation of F*(s) can be simplified by a substitution of variable (1, 3,4) 
i= rial which gives, 

F*(s) = F(z) (5) 
This change of variable maps one strip of width jWs on the s-plane into the whole 
z-plane. Consequently, F(z) has the same number of poles as F(s) and approxi- 
mately the same number of zeros. F(s)can now be evaluated by either a geometric 
or a potential analog of the z-plane and, if desired, the results can be transformed 
back into the s-plane or the time domain. This inverse transformation is usually 
not necessary because the usual criteria of stability, damping ratio, minimum set- 
tling time, and frequency response can be determined directly from the z-plane. 

There are certain advantages in using the log z-plane instead of the z-plane 
with a potential-plane factor analog (2). 

Some of these advantages are as follows: 

(1) Several decades of the variable may be used. 

(2) The percentage accuracy is constant over the whole log z-plane. 

(3) Edge effects are easily minimized. 

(4) Two different types of input information may be used, (z = Zi) or be =) ‘ 

(5) Scale expansion is easily obtained by a change of variable, e.g. 

(z'=z? or z'=z-2Z)) 
(6) Most performance criteria are straight lines on the log z-plane. (Damping 


ratio, minimum settling time, etc., are straight lines.) 
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FIGURE 3 A COMPARISON OF (a) THE s-PLANE, (b) THE z-PLANE, AND 
(c) THE LOG z-PLANE 


Figure 3 shows the relationships between the s-plane, z-plane, and log z-plane. 
The left half of the s-plane maps into the area inside the unit circle in the z-plane 
and into the area to the left of the Izl=I line in the log z-plane. A line of constant 
damping ratio ¢ in the s-plane maps into a heart-shaped curve in the z-plane and 
into a straight line in the log z-plane./ 

A derivation of the simple relationship between the s-plane and the log z-plane 
is shown below: 


sT _ site! _ (oe +ju)T 
oe = 


z=e (6) 


Log. (z) = sT = IsiTe/45 = oTt+jwt (7) 


where T is the sampling interval which is a positive real number. Equation (7) 
shows that rectangular coordinates on the log z-plane are the same as the rectangu- 
lar coordinates on the s-plane with the scales multiplied by T. Consequently, the 
s-plane and log z-plane plots will differ in appearance only by a scale factor T. 


Caution must be exercized with respect to the poles or zeros at $ =—@ since 
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these poles or zeros will appear at the origin of the log z-plane. These poles or 
zeros must be taken into account when evaluating F(z). 
The Pulse Transfer Function: 

In linear system analysis, the transfer function concept is very useful for 
relating the output response of the system to the input excitation to the system. For 
continuous linear systems, the transfer function is the Laplace transform of the out- 
put response with an impulse input excitation. In linear sampled-data systems, 
there is an analogous concept which is called a pulse transfer function. The pulse 
transfer function relates the output pulse sequence to the input pulse sequence and 
is valid only at the sampling instants. The pulse transfer function gives no indica- 
tion of the behavior of the system between sampling instants. 

The Z-Transform: 

In linear system problems, a continuous function of time f(t) can be trans- 
formed into an algebraic function of a complex variable F(s) by using the Laplace 
transform concept. In linear sampled-data problems, a discrete function of time 
f*(t) can be transformed into a function of a complex variable F(z) by using an 
analogous concept which is the Z-transform. There is also a relation between 
F(s)and F(z) where z =e! , such that F(z) can be found from a given F(s)and 
vice versa. Tables of Z-transforms are available in the literature (3,4) which 
usually show three columns, f(t) , F(s), ana F(z). 


For Example: 


eT 
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The procedure for finding the Z-transform F(z) froma given F(s) is 
outlined below: 


(a) Expand the F(s) by a partial fraction expansion to get 


SP | ec Lire eae 
F(s) = K AS a (s—s2) * (e=s3 + (8) 


(b) Evaluate the residues A, B,C, °** , (these can be found with a potential- 
plane factor analog, or often by inspection). 


, S 
(c) Using the same residues, convert the s-poles to z-poles by letting zj =e 


=p 
and set 
To. = gem = 9 pean ee 
eaeayeD ; 
Z Zz ay & 


(d) Cross-multiply to find a common denominator which will be the product 
of the z-poles. 


(e) Factor the numerator to get the z-zeros. (A potential-plane analog can 
be used for this operation). 


The result will be a function F(z) in the form 


F(z) = K Res 


If an ESIAC® factor analog computer (2) is used, factors of the form (z = Zi) or 


| Zi 
( — jar more desirable than the form ( — 3) shown in Equation (10). The 
eo 


factor type (z a zi) can be obtained from Equation (10) by multiplying by 


z. 
I nN 
Zz 
Np 
Z 
ar 
Pp gP 


®*Trademark registered by Electro Scientific Industries, Inc., Portland, Oregon 
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to get 
Wap (2-22) 
F(z) = K a PT (11) 
HI Np (z—2,) : 


The effect of the change in factor type results in a change in the order of the 
pole or zero at the origin. If all of the z-poles are first order, the result of the 
factor change adds a zero at the origin. This can be easily seen by factoring out 
a zin Equation (9). 

An example of the use of the Z-transform is shown in the derivation of G(z) 
in the Appendix. 

ILLUSTRATIVE EXAMPLE: 

A simple example of a sampled-data control system is shown in Figure l. In 
the following discussion, this system is analyzed with an ESIAC potential-plane 
analog to show the effect of adding a compensating network. The system response 
without the compensating network is found first as shown below. 

The pulse transfer function G(z) is found from the corresponding G(s) by using 
the theory of the Z-transform as shown in the Appendix. 

From the Appendix: 


0.368K (z+ 0.72) 


=" (z—1) (z-—0.368) 


(12) 
Then the closed-loop response is 


I tg OORT RT) 
R(z) 1+G(z) — (z—1) (z— 0.368) + 0.368K (z +0.72) (13) 


Clz 
The poles of e are found by factoring the denominator of Equation (13). 
Z 
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First set the denominator equal to zero and rearrange as follows: 


a) 


(7-1) (z-—0.368) + 0.368K (z+0.72) = 0 


or 


in 
= 
| 
ll 
aw | ~ 
“8 
oO} mM 
+1o 
Sy i 
eee 
ee) abs 
oO;—= 
as I 
MIN 
Or 


This equation is easily solved by plotting a root-locus plot on the log z-plane 


with an ESIAC potential-plane analog as shown in Figure 4. 
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FIGURE 4 ROOT LOCUS PLOT ON THE LOG-z-PLANE. 
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SAMPLED DATA SYSTEM 


WITHOUT COMPENSATION 


The unstable region where Iz| >| corresponds to the area outside of the unit 


For a damping ratio of 0.5, K=0.57 and the roots of 


circle on the z-plane. 


Equation (15) are 


z= 0.72236°, 0.721-36° 


The maximum value of Kp is found at the point where the loci crosses into the 


unstable region. At this value, Km = 2.40, the system will begin to oscillate. 


This system may be improved by adding a compensating network which, in 


For Example: 


effect, changes the root locus curve to some more desirable shape. 


it may be desirable to obtain more gain for the same value of damping ratio. 
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(18) 
(19) 


Let us adda zero at 
For a damp- 
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The F(z) shown in Figure 1 is then 
(z-0.1) 
F(z) G(z) 

1+ F(z) G(z) 


LEOe 
(0.2) (z—0.5) 


F(z) = 
R(z) 
F(z)G(z) = -l 


C(2) 
(z—0.1)(z —0.368)(z—1) 
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The effect of adding a pulsed lead network to the system can be shown by adding 
K(0.2) (0.368) (z—0.5) (z+ 0.72) 


A physical network that gives this F(z) can be found by using the method shown 


The closed-loop response with the pulsed lead network added is 


The root locus plot of the above equation is shown in Figure 5. 
ing ratio of 0.5,K=5.2 compared toK=O,5/ for the uncompensated system. This 


a pole-zero pair to Figure 4 and replotting the root loci. 
z = 0.5£0° anda pole atz = 0.1Z0°. 


in Reference 3, pages 136-43. 
and the characteristic equation is 


or 
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FIGURE 5 ROOT LOCUS PLOT ON THE LOG-z-PLANE. 
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increase in gain corresponds to an improvement in system performance since 
the system error is inversely proportional to open loop gain. 
The closed-loop pole locations are found from Figure 5 and are at 


z= 0.54Z61°, 0.542-61°, 0.6120° 


The closed-loop transfer function is 


R(z)  (z-0.54Z261°) (z-0.54 Z—61°) (z— 0.61) oy) 
C(z) 


Curves of the magnitude and angle of versus log frequency can be plotted 


R(z 
if desired. This corresponds to the Bode plots for continuous systems in which 


| 2 
$= jw. The data for the log frequency plots is obtained by evaluating —— for 


R(z) 
lzl=l and Zz varied from 0 to 180°. The relationship between real frequency f and 
fz is 
emf 
Zz=wT= ; (21) 
s 


The data for the log frequency plots will be periodic with respect to one-half 
of the sampling frequency. This illustrates the ambiguity that results when the 
sampling frequency is not at least twice the highest frequency of interest in the 
control system. 

SUMMARY: 

Some fundamental concepts of sampled-data systems are presented, and 
the relationships between the s-plane, the z-plane, and the log z-plane are 
shown. The advantages of the log z-plane over the z-plane are demonstrated 
by an illustrative example calculated with the aid of a potential-plane factor 


analog computer. 
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APPENDIX 


EXAMPLE OF Z-TRANSFORM CALCULATION 


= a K 

Given: G(s) —— es ea ene 

S s(s+a) 
Required: 
Find the Z-transform of G(s). 
First expand G(s) into a partial fraction expansion and get 
( am) A B C 
l—e = 
S sta 


G(s) = 


The coefficients A, B, C, can be evaluated with the ESIAC”, but for this 


Simple example they can be found by inspection. 


2. 2 


ie) K K K 
as a°s  —a*(s ta) 


G(s) = 


Now take the Z-transform of each term individually and get: (see Reference 3 


or 4 for Z-transform table) 


G(z) = Gea K lt, Kt. 4H 2 
a (z-1) a” (z-lI) ao 2z-e 
where Z = ean eee = 271 
z Z 
Then G(z) = x ete + —_i-— 
a} (z—1) a a(z—e ) 


Now find a common denominator for the term in brackets and get: 


K_| Ta eae) + Cm Caen 


° a(z—I) (2— aap 


For this problem, T=a=I, and G(z) reduces to 


G(z) = 


K (2+0.72) 0.368K (z +0.72) 
Gz) = amy ety 
(z-i)(z—e ) (2-)(z-e ) 
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